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Abstract 

We investigate the thermodynamics of a Bose gas interacting with repul- 
sive forces and confined in a harmonic trap. We show that the relevant pa- 
rameters of the system (temperature, number N of atoms, harmonic oscillator 
length, deformation of the trap, s-wave scattering length) fix its large N ther- 
modynamic behaviour through two dimensionless scaling parameters. These 
are the reduced temperature t = T/T^ and the ratio rj between the T = 
value of the chemical potential, evaluated in the Thomas-Fermi limit, and the 
critical temperature T^ of the non-interacting model. The scaling functions 
relative to the condensate fraction, energy, chemical potential and moment of 

inertia are calculated within the Popov approximation. 
02.70.Lq, 67.40.Db 
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First measurements of relevant thermodynamic properties in gases of alkali atoms con- 
fined in magnetic traps have been recently become available . These include the temper- 
ature dependence of the fraction of atoms in the condensate and the release energy. On the 
theoretical side the study of non uniform interacting Bose gases has a long history, starting 
from the works by Gross, Pitaevskii and Fetter |^ and the first approaches to thermody- 
namics developed in ref. Q]. More recently several theoretical papers have focused on the 
problem of the thermodynamic behaviour of such systems fsHlOl, and first calculations of 



the temperature dependence of the condensate have become available 

The properties of these trapped Bose gases at a given temperature are fixed by a large 
number of parameters: total number of atoms in the trap, atomic mass, harmonic oscillator 
length, deformation of the trap and strength of the interaction, which in a dilute gas is fixed 
by the s-wave scattering length. This makes these systems very rich for both experimental 
and theoretical investigations, but at the same time difficult to study in a systematic way. 
For example the deformation of the trap in the experiments carried out at MIT and at 
Jila as well as the corresponding number of atoms are quite different and it is not always 
easy to compare the experimental results. The purpose of this letter is to point out the 
occurrence of a scaling behaviour in the thermodynamics of these many-body systems. This 
permits to compare very different experimental situations, provided they correspond to the 
same value of the scaling parameters. 

A first natural scaling parameter is given by the reduced temperature t = T/Tj^, where 
kj^T^ = huj{N/CiS)y/''^ ini, with C(3) ~ 1.202, is the critical temperature predicted by 
the noninteracting harmonic oscillator model in the large N limit and u = (cOxtOyCUzY^^ is 
the geometrical average of the oscillator frequencies characterizing the confining potential 

Vext = \m{ujlx'^ + ^2y2 ^ ^2^2) _ 

In addition to the reduced temperature t we introduce a second scaling parameter ac- 
counting for the role of the two-body repulsive interaction. This parameter is fixed by the 
ratio 7] between the T = value of the chemical potential, calculated in the Thomas-Fermi 
approximation, and the critical temperature fcsT^. In the Thomas-Fermi approximation 



the ground state density of a Bose gas takes the form [0 no(r) = (/Uq^ — Vextij))/ g, where 
g = Anh^a/m is the couphng constant fixed by the (positive) s-wave scattering length a and 
the chemical potential is given by /i^^ = ^^^ujilhN a/ auoY^^ ^ where auo = {h/mcuY^'^ is 
the harmonic oscillator length. The parameter rj then takes the useful form 




with a = ((3)^/^15^/72 ^ 1-57. Note that r] exhibits a very smooth dependence on N 
(rj ~ N^^^^) and its value in the presently available experiments ranges from 0.38 to 0.45 
even if N varies by orders of magnitute in the different experiments. The ratio (0) can be 
also written in the form r] = 2.24(a^no(i" = 0))^/^, where no(r = 0) is the density of the 
atomic cloud in the center of the trap evaluated at zero temperature. With the above values 
of rj, the gas parameter a^no{r = 0) turns out to be extremely small (10~^-10~^) so that the 
system is always very dilute. 

In the following we will show that the relevant thermodynamic properties of the system 
exhibit, when expressed in their natural units, a scaling behaviour in terms of the parameters 
t and 77. Scaling is actually achieved only if N is sufficiently large and finite size corrections 
to the asymptotic behaviour of the various thermodynamic functions can be ignored. We will 
show that scaling is reached with very good accuracy for most thermodynamic properties 
not only in the experimental conditions of ref. where N is of order of 10^ — 10'', but also 
in the ones of ref. where is much smaller (10^ — 10^). 

The scaling behaviour can be exploited by investigating the thermodynamic behaviour 
of these systems within the Popov approximation |l^J^,p]. This approximation provides 
coupled equations for the condensate density and the thermal density of the system. In 
this letter we will use the semiclassical approximation for the excited states and hence we 
will always consider temperatures much larger than the oscillator temperature huj/ks- The 
equation for the condensate is given by 

- — V2$(r) + [Vextir) + g{no{r) + 2nT(r))]$(r) = 0, (2) 
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and the dispersion law of the elementary excitations takes the Bogoliubov-type form 
e(P,r) 



^ + K.i(r) - /i + 2gniv)^ - g^nliv). (3) 



In eqs. (@), (H) rio = is the density of the condensate, n = + no is the total density 
of the system and the thermal density ut is defined by nT{r) = J dp/ {2nh)'^ {-de/ dfi) f (Pe) 
where /(/3e) = {exp{(3e) — 1)^^ is the quasi-particle distribution function. The semiclassical 
approximation (H) to the excitation spectrum is accurate only if the non-condensate density 
riT varies smoothly on the scale of the oscillator length auo-, ^^nd if the relevant values of 
p satisfy the condition p ^ Tijauo-, which also implies e ^ Tiuj. The above approximation, 
called Popov approximation, coincides with the Bogoliubov approach at very low temper- 
atures, while at high T it coincicides with the finite temperature Hartree-Fock theory for 
Bose systems (for further details on the Popov approximation see refs. [§],§). Equations 
(@)"(i) si^re solved using a self consistent procedure, the value of the chemical potential \x 
being determined by the normalization condition / (ir(rio(r) + ^^(r)) = A^q + = ^• 

The scaling behaviour of the Popov equation is exhibited in the Thomas- Fermi limit, 
which is reached when the ratio fs^^ /fiuj = 0.94:N^^^r] is much larger than unity. In this 
limit, which is formally achieved taking N ^ oo with t] fixed, the kinetic energy term in 
the Schrodinger equation (0) for the condensate can be neglected. By introducing the di- 



mensionless quantities pi = ^Jl/2mkBT^p^, fi = ^m/2kBT^uJiri, ji = fi/ksTl^, e = e/kBTj^, 
fiT = nT{2kBT^ / moj'^Y/'^ / N and no = no(2kBT^ / mcu'^)^/'^ / N , one can rewrite the Popov 
equations in the reduced form 



i (/i - f2 - 2ghT) 9{fi - f2 _ 2ghT) , (4) 



no(r) 



e(p, r) = V (p2 + r2 - /i + 2g{no + ut)) - g'^nl , (5) 

where hrif) = l/(7r3C(3)) / rfp(-9e/9/i)/(e/t), g = Mrj^l'^/lh = h.2N^/^a/aHo and e{x) is 
the step function. The normalization condition for the reduced densities reads: / (if (no + 
fiT) = 1- Equations (§)-(0) exhibit the anticipated scaling behaviour in the variables t and 



rj. Starting from their solutions one can calculate the condensate and the thermal densities 
as well as the excitation spectrum. This gives access to all the relevant thermodynamic 
quantities of the system. For example the condensate fraction is given by 

1 - ^ / drdp (-^] fm . (6) 



N 7r3C(3) J \ dfL^ 

In particular the vanishing of the right hand side of eq. (j^) fixes the value of the critical 
temperature tc = Tc/T^. To the lowest order in the renormalized coupling constant g 
one finds ||8|JT0[| t^. — 1 — 0.437]^/^ = 1 — \.?)N^^^a/ ano- The energy of the system can be 
calculated starting from the standard thermodynamic relation in terms of the entropy. In 
units of NksT^ one finds: 

E 5 r , ,ds 



where the entropy per particle s is given by the combinatorial expression 

sit, v) = / dHp[{i + fm) iog(i + fm) 

-fm^ogfm] ■ (8) 

Another important quantity is the normal (non superfluid) density given by the natural 
generalization 



of the Landau formula |TB| to the case of non uniform systems. Note that in general the 
normal density p„ differs from the non condensate density rij-. Only at temperatures larger 
than the chemical potential one has p„ ~ mn-r- In terms of p„ one can calculate the moment 
of inertia, defined as the linear response function to a crancking rotational field [jl6|. In the 
case of an axially symmetric trap one finds the result 

9 = 1 dr{x^ + y^)pn{r) . (10) 

Deviations of G from the rigid value Qrigid = I c?r(x^ + y'^)n{r) provide a signature of the 
superfluid behaviour of the system. Above one has p„ = mriT = mn and O = Qrigid- 



In Fig.l we explicitly show the accuracy of the scaling behaviour, by plotting the scaling 
function (^) calculated with rj = 0.45 and the results obtained from the full solution of the 
Popov equations with different choices of and of the other parameters, yielding 

the same value r]. The open circles have been obtained with the choice: a/ano = 7.35 x 10^^, 
A = uJz/uJx = V8, N = 5 X 10^, which corresponds to the experimental conditions of ref. 0, 
while the solid circles have been calculated with the choice: a/ano = 2.55x 10^^, A = 18/320, 

= 2.9 X 10^, which is close to the experimental situation of ref. This example shows how 
extremely different experimental conditions give rise to the same thermodynamic behaviour. 
The figure clearly shows that scaling is very well verified for these configurations. Only very 
close to Tc the = 5 x 10^ points exhibit small deviations from the scaling behaviour. In 
fact close to the scaling law (|^) is approached more slowly with increasing A^. In the 
same figure we also plot the predictions for the same configurations obtained switching off 
the two-body interaction. Deviations from the non-interacting result 1 — are due to finite 
size effects, whose importance, in the presence of the interaction, turns out to be strongly 
quenched. The scaling behaviour has been verified also for the other values of the parameter 
1] and for all the other thermodynamic properties considered in this work. 

In Fig. 2 we present results for the condensate fraction Nq/N as a function of the reduced 
temperature t for three different values of the scaling parameter rj, covering the presently 
available experimental conditions. The open diamonds with the error bars are the results 
of the Monte-Carlo simulation of ref. which correspond to the value 1] = 0.33 and which 
are in good agreement with our predictions. The dots are the experimental results of ref. 
0. In the experiments the number of particles A^ varies with T, with the value of rj ranging 
from 0.45 to 0.39. The experiments exhibit smaller deviations from the non-interacting 
curve with respect to the theoretical predictions. One should however keep in mind that the 
measured value of T corresponds to the temperature of the thermal cloud after expansion. 
The identification of this temperature with the one of the system before expansion ignores 
the interaction with the condensate which is expected to produce an acceleration of the 
thermal cloud. A preliminary estimation shows that for T ~ 0.5T° the final kinetic energy 
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of the thermal cloud is about 10% larger than its value before the expansion. 

In Fig. 3 we present the results for the chemical potential in units of ksT^ corresponding 
to the same values of 1]. Notice that for t —>■ the plotted quantity coincides with rj (see 
eq. (|1|)). In the classical limit, T ^ T°, the dependence on the interaction parameter rj 
disappears and one finds the classical ideal gas prediction n/ksT'^ = t log(C(3)/t'^). 

In Fig. 4 we report the results for the total energy E. At high temperature the behaviour 
is given by the classical law E/ksT^ = 3t. In Fig. 4 (inset) we also show our predictions 
for the release energy E^, defined as the energy of the system after the trap potential has 
been switched off. In the classical non-interacting limit one has Efj = 3NkBT/2. First 
experimental results for the release energy have been obtained in ref. and are shown in 
the figure. 

Finally in Fig. 5 we show our predictions for the moment of inertia. We find that, dif- 
ferently from the other quantities discussed above, the ratio Q/Qrig does not exhibit any 
significant dependence on rj. In the same figure we also show the moment of inertia calcu- 
lated in the non interacting model for two different values of and deformation A. It is 
worth noting that in the absence of interaction the value of G/ Qrig depends rather crucially 
on N and A. 

In conclusion we have explored the scaling behaviour of the thermodynamic properties of 
a trapped Bose gas. The new key parameter which permits to discuss the role of two-body 
interactions is the scaling parameter (|l]), given by the value of the zero-temperature chemical 
potential in units of ksT^. Scaling is predicted to work very well in the typical conditions 
of available experiments. 

We would like to thank D. S. Jin for providing the experimental results on the condensate 
fraction and release energy 0. One of us (S.G.) would like to thank the hospitality of JILA 
while part of this work was performed. 
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FIGURES 

FIG. 1. Condensate fraction as a function of T/T^ for r] = 0.45 (solid line). The open circles 
refer to = 5 x 10^ Rb atoms in the JILA-type trap. The solid circles correspond to = 2.9 x 10^ 
Na atoms in the MIT-type trap. The dotted line is the 1 — curve of the non-interacting model 
in the large limit. The open and solid triangles correspond to = 5 x 10^ and A^ = 2.9 x 10^ 
non-interacting particles in the JILA and MIT-type traps respectively. 

FIG. 2. Condensate fraction as a function of T/T^ for three values of the scaling parameter 
r]. Solid line: r] = 0.45, long-dashed line: r] = 0.39, short-dashed line: r] = 0.31. Open diamonds: 
PIMC results of ref. Solid circles: experimental results from ref. Q. The dotted line refers to 
the non-interacting model in the large A^ limit. 

FIG. 3. Chemical potential as a function of T/T^ for three values of the scaling parameter r] 
(see Fig. 2). The dotted line refers to the non-interacting model in the large A^ limit. 

FIG. 4. Total energy of the system as a function of T/T^ for three values of the scaling pa- 
rameter r] (see Fig. 2). (inset) Release energy for the same values of ry. The solid circles are the 
experimental results of ref. Q. 

The dotted lines refer to the non-interacting model in the large N limit. 

FIG. 5. Ratio @/&rigid as a function of T/T^ for three values of the scaling parameter rj (see 
Fig. 2). The three curves coincide almost exactly and are represented by the solid line. The 
dotted line refers to A^ = 5 x 10^ atoms in the JILA-type trap in the non-interacting model. The 
dot-dashed line refers to A^ = 2.9 x 10'' atoms in the MIT-type trap in the non-interacting model. 
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